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Abstract 

We present the consequences of a large radiative correction term coming from 
Supersymmetry (SUSY) upon the electron neutrino fluxes streaming off a core- 
collapse supernova using a 3-flavour neutrino-neutrino interaction code. We ex- 
plore the interplay between the neutrino-neutrino interaction and the effects of 
the resonance associated with the fj, — r neutrino index of refraction. We find 
that sizeable effects may be visible in the flux on Earth and, consequently, on the 
number of events upon the energy signal of electron neutrinos in a liquid argon 
detector. Such effects could lead to a probe for Beyond Standard Model (BSM) 
physics and, ideally, to constraints in the SUSY parameter space. 
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1 Introduction 



During the past decade, our comprehension of the neutrino interactions in a supernova 
environment has deeply evolved, due to the increasing precision of neutrino experi- 
ments and computational force. In addition with the traditional MSW effect [U 12], the 
neutrino-neutrino interaction has been proven to be fundamental for the neutrino evo- 
lution in such environment. This interaction, whose effects have changed the existing 
paradigm of SN neutrino physics, has been recently the subject of intense investigation 



[31 SI El El [71 El El [101 [HI H2l H3l HH H5l H6l HTl HH H9l EDI [SI]. Fora recent extensive 



review on this subject, see [22] . 

In addition with these interactions, an other matter interaction V^ T , which arises at 
one- loop order (O ( — ■ "a ttz~ )) makes the mu and tau neutrino index of refraction 



different (231 121] • This radiative correction potential is related to the charged- current 
matter potential V c , in an electrically neutral medium, by a parameter e, which is 
defined as —p u An Tfl , = V^ T = eV^i The value of e is 5.4 x 10" 5 in the Standard Model. 
Only a few papers have shown the interest of such radiative correction term. Without 
the neutrino-neutrino interaction, V^ T has been shown to present effects on Earth if the 
Up and u T fluxes are different at emission [25]. Via a non-zero CP- violating phase, V^ T 
can also in theory induce effects on the (anti-) electron neutrino fluxes in the Sun [26] or 
in a supernova environment [271 12H] • Including the collective effects, the CP- violation 
effects due to V^ r turn out to be larger [29], and the electron neutrino flux displays a 
#23 dependence in presence of a large V^ T coming from a high density profile [301 EI] • 

Despite its success, the Standard Model (SM) of particle physics exhibits a number 
of shortcomings which might be remedied by new Physics showing up at the Terascale. 
One very popular extension of the SM is Supersymmetry (SUSY) [32]. Virtues of SUSY 
are numerous: it can provide a natural candidate to explain dark matter and it allows 
the elimination of the hierarchy problem and the unification of the gauge couplings at 
the scale of grand unification (M GUT ). 

In a previous paper, we have calculated all one-loop radiative correction terms to 
the matter interaction in the SUSY framework [33] . SUSY particles like sleptons, 
squarks up and down, charged Higgs, neutralinos and charginos take place in these 
loops. With a dedicated numerical routine, we have scanned and identified regions 
in the parameter space that yield interesting value of the V^ T potential, namely large 
values up to Ssusy — 2.10~ 2 . 

The recent simulations of neutrino evolution [311 [35] in supernova take into account 
a dynamical density profile, in addition with the neutrino collective effects. All in- 
teractions combined show a characteristic imprint in the anti-neutrino fluxes, whose 
detection could yield precious information concerning the dynamics of the density pro- 
file or fundamental neutrinos properties like the hierarchy or the value of the third 
mixing angle 6\%. With larger value of e than in the SM case, we show that V^ T can 
exhibit sizeable and characteristic effects on the electron neutrino fluxes particularly 

1 We do not consider the presence of other charged leptons than electrons in the supernova envi- 
ronment. Therefore, V c = V e . 
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in the inverse hierarchy where the v — v interaction display salient features such as 
the synchronization regime, the bipolar transition and the spectral split. Such BSM 
imprints could be seen in a neutrino observatory on Earth and could ideally lead to 
constraints in the SUSY parameter space. 

The paper is organized as follows. In Sec. [2] we explicit the evolution equations 
which lead to the calculation of the electron neutrino flux, describe the framework 
used, and justify the approximations taken to yield our results. In Sec. [3] we describe 
the three typical behaviours that can occur due to the interplay between the neutrino- 
neutrino interaction and the \i — r resonance. In Sec. H] and Sec. [5j we analyze the 
role of the parameters that can influence such interplay and consequently the v e flux 
as a function of energy. The subsections 15.21 and 15.31 are dedicated to the study of the 
consequence of varying the luminosity and the density which is equivalent to looking at 
the flux at different times. We then use these fluxes to display the number of events in 
a liquid argon detector that we describe in Sec. [6j Finally, before concluding in Sec. [H 
we discuss in Sec. [3, about the possibility of the observation of such beyond standard 
effects for realistic conditions. 

2 Theoretical framework 

The main goal of this paper is to explore the impacts of the parameters that influence 
the interplay between the neutrino-neutrino interaction and the \i — r resonance and 
consequently study the possibility that such effects can be observable in a neutrino 
observatory in order to get a possible probe for models BSM. In this section, we first 
present the evolution equation we have numerically solved and the hypothesis used to 
yield the results of the following sections. 

2.1 Equation of evolution. 

In a dense environment the non-linear coupled neutrino evolution equations with neu- 
trino self-interactions are given b}0: 



where ip Va denotes a neutrino created at the neutrinosphere initially in a flavour state 
a = e,fi,r, H = UH vac U^ is the Hamiltonian describing the vacuum oscillations 
H vac = diag(£' 1 , E 2 , E 3 ), E i= i^ >3 being the energies of the neutrino mass eigenstates, 
and U the unitary Maki-Nakagawa-Sakata-Pontecorvo matrix 



/I \ / C13 s 13 e~ iS \ / c 12 s 12 0\ 
U = r 23 T 13 T 12 = C 23 S 23 1 -Sia c 12 , (2) 

\0 -S 23 C23/ \S 13 e i5 C13 / V 1/ 



2 The dependence on r is equivalent to the dependence in t since we suppose that neutrinos travel 
at light speed in the supernova and we take c = 1 . 




(1) 



2 



C{j = cosOij (sij = sin%) with #12, #23 and #13 the three neutrino mixing angles. The 
presence of a Dirac 5 phase in Eq.(J2]) renders U complex and introduces a difference 
between matter and anti-matter. 

The neutrino interaction with matter is taken into account through an effective 
Hamiltonian which corresponds to the diagonal matrix H m = diag(V e , 0, V^ T ). V e (r) = 
V2Gpn e (r) is the matter potential due to the charged- current interaction between 
electron (anti-) neutrinos and the electrons present in the medium where Gp is the 
Fermi coupling constant and n e (r) is the electron density in the star. At the tree level, 
neutral current interactions introduce an overall phase only. Note that this density will 
imply in supernova possibly two resonances when the MSW condition [2] (here shown 
for two flavours) is reached: 

— cos26 = V e (r) (3) 

Depending on the sign of Am| 2 ,the " H-resonance" (high density) may happen either for 
neutrinos or anti-neutrinos while the " L-resonance" (low density) happens for neutrinos 
only and is dependent on Am^. 

As previously mentioned, in the Standard Model case, Botella et al. in [23] have 
highlighted the presence of a one-loop matter potential arising from radiative correc- 
tions to neutral-current and v T scattering. Such matter potential can therefore be 
seen as an effective presence of r particles and one can write 

V llT = V2G F Y T cS SM n B (4) 

where n B is the baryon density inside the supernova and 



s 3V2G F m 2 T 

X t,SM - 



(27T) 
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2.7 xl(T 5 , (5) 



The value of Yfg M is obtained using the hypothesis of a isoscalar medium i.e, n e = 
n p = n n {rip and n n being respectively the proton and neutron density in the star), 
therefore the neutron fraction Y n = n n / {n p + n n ) is equal to 0.5. For our convenience, 
we have defined the overall radiative correction factor as : 

VeS _i_ vcff 

I t,SM ~r I t,SUSY /„n 

£ = y (6) 

1 e 

where F r e f USY = in the Standard Model and the electron fraction Y e = n e / {n n + n p ) 
is taken to be 0.5 all along the paper. 

The general form of the neutrino self-interaction term is 



(7) 



where p = p Ua (—p* Ua ) is the density matrix for neutrinos (antineutrinos) 

P^= $tfl>»* I^J 2 
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q (q') denotes the momentum of the neutrino of interest (background neutrino) and 
dn a is the differential number density. The emission geometry is based on the so 
called "bulb model" with spherical symmetry [12] : neutrinos are assumed to be half- 
isotropically emitted from the neutrinosphere. 

In the multi-angle treatment, SN neutrinos travel on different trajectories. This is 
taken into account by the factor (1 — q • q') in Eq.([7|). The single angle approximation 
consists of averaging such factor along the polar axis, i.e. p(q) = p(q), Eq.© reduces 
to 



2-kRI 

with the geometrical factor 



(9) 



D(r/R) = - 



(10) 



where the radius of the neutrino sphere is R = 10km. In a spherically symmetric 
environment, a single- angle treatment tends to capture the main non- linear behaviour. 

For this paper, we shall adopt 

Ef, 



Tl (E U JF 2 ( V ) 1 + exp (E„JT^ - V ) 



as the neutrino luminosity, where -^2(77) is the Fermi integral, L° and T Va are the 
luminosity and temperature at the neutrinosphere. (E Ua ) is the average neutrino energy 
of the corresponding flavour a. 

As we are going to see in the following sections, the neutrino-neutrino interaction is 
characterized by three distinct phases for the neutrino evolution: the synchronization 
regime, the bipolar transition and the spectral split. The synchronized regime takes 
place in the first 50 km outside the neutrino sphere. In this phase, the strong neutrino- 
neutrino interaction makes neutrinos of all energies oscillate with the same frequency 
so that flavour conversion is frozen [HI [15]. When the neutrino self- interaction term 
diminishes, large bipolar oscillations appear that produce strong flavour conversion for 
both neutrinos and anti-neutrinos for the case of inverted hierarchy, as long as the 
6*i3 value [13] is non-zero. Eventually neutrinos show complete (no) flavour conversion 
for energies larger (smaller) than a characteristic energy E c ~ 8MeV which can be 
related to the lepton number conservation [17]. This is the spectral split phenomenon: 
while electron neutrinos swap their spectra with muon and tau neutrinos; the electron 
anti-neutrinos show a complete spectral swapping. Such behaviours are found for both 
large and small values of the third neutrino mixing angle, contrary to the standard 
MSW effect. 

We have chosen here to describe the evolution equation using wave-functions such 
as in [12j [29] for instance. Note that equivalently the density matrix or the polariza- 
tion vector formalism can be used to describe the neutrino evolution in the supernova 
environment. 
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2.2 Work hypothesis. 



We give the general framework and the hypothesis that we make to obtain our numerical 
results. 

The results we obtained use the best fit oscillation parameters to date, i.e. Am^ = 
7.6 x 10" 5 eV 2 , sin 2 2#i2 = 0.87 and | Am 2 ^ | = | Am 2 atm \= 2.43 x 10~ 3 eV 2 for the solar 
and atmospheric differences of the mass squared! and mixings, respectively [36] • Apart 
from Sec. 15. 1\ we choose #23 = 40° in the first octant, i.e equivalent to sin 2 2^23 = 0.97. 
Unless notified, the Dirac phase is taken to be zero and sin 2 2^13 = 0.1. This value of 
#13 taken is in agreement with the current experimental upper-limit [36] . 

During a supernova explosion, three distinct phases take place: the prompt delep- 
tonization burst, the accretion phase, and the cooling phase. In this paper, we will 
focus on the neutrino signal produced during the early cooling phase where the density 
of neutrinos is sufficiently important to induce sizeable non-linear behaviour. Concern- 
ing the cooling phase, the hypothesis of equal luminosities and equipartition of energies 
among all neutrino flavours at the neutrino-sphere is made for the neutrino signal upon 
which we are focusing. We also make the assumption of an exponential decrease of the 
luminosity 

L v = L vfi x exp(-t/r) (12) 

with L VQ = 10 52 erg • s _1 and r = 3.5 s the typical cooling time. Apart from Sec. 15.21 
and Sec. |6j we look at time t ~ 8s after post-bounce which corresponds to L u = 10 51 
erg • s -1 . We consider the hierarchy i.e. (E Ve ) < (Ep e ) < (E Ux ) with typical values 
of 10, 15 and 24 MeV respectively. These values correspond approximatively to the 
cooling case I in [37]. 

In this paper, one of the observables is the v e flux, which can be written 

= P{v e ^Ve){F*-F«) + F« (13) 

using the unitarity of the evolution operator of Eq.([T]) and where F^ is the flux of 
flavour a emitted initially. We here made the hypothesis of equal and v T initial 
fluxes. We can link it to the luminosity L Ua (E u ) defined in Eq. (jTTj) via the relation: 

Moreover, we define the flux received on Earth by 

^Earth =< ^ > = < p[Ve ^ ^ >cMt {p R _ pRj + p R ^ 

The probability exiting the SN < P[y e — > u e ) > cx a, and therefore the flux on Earth, 
is averaged over distance because of the decoherence phenomenon that takes place 
for neutrinos travelling over large distance in vacuum [38]. As the neutrinos stream 



2 2 A 2 

5 Note that E t - E j 
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off the supernova, a shock wave propagates through the SN envelope coming with a 
reverse shock. This induces a dynamical density profile which leads to the breaking of 
the adiabaticity of the H-resonance and possibly to a multiple resonance effect. Such 
interplay with the neutrino fluxes may imply a salient fingerprint on the signal received 
in an observatory on Earth as pointed out in [521 EH SD EH] • 

All along the paper we will focus on the electron neutrino flux in the inverted hier- 
archy. Consequently, the shock-wave effect should not impact on the electron neutrino 
flux since no H-resonance occurs. Therefore, for our electron number density, we adopt 
the usual simple analytical profile for the : 

A(r) = V2G F n B = X r f^j . (16) 

where R is equal to 10km. Apart from Sec. 15.31 and Sec. [61 we take for reference 
A r = 5.7 x 10 5 km~ 1 which corresponds to a late time density profile. 

In [12], SM one- loop corrections to v — v interaction have been studied and turned 
out to be of the same order as Y* SM or even smaller. Though they have not been 
calculated yetQ, supersymmetric one- loop corrections iov—v interaction will be similar 
to the SM ones but only larger in the same way than for USY in comparison with 
Yt%m- For £ < 2.10" 2 , we expect them to be negligible. All the hypothesis made in 
this section will be discussed in Sec. [3 



3 Typical behaviour of a large fx-r radiative correc- 
tion 

In this section we review the consequences of a large V^ T on the v e flux in the inverted 
hierarchy, in presence of the neutrino-neutrino interaction. Unlike in [30], we do not 
increase the value of the density A r to obtain a large effective r matter potential but 
we use instead the radiative correction term e whose value can rise depending on the 
SUSY parameters [33] . Consequently, varying the e for a given energy will only move 
the fi — t resonance above the synchronized region leaving unchanged the position of 
the L-resonance. 

In Fig.([T]) we show the oscillation probabilities for three different energies E v = 15, 
30 and 35 MeV with e = e r = 10~ 3 in order to identify precisely the physical effects 
appearing. Analogously to the H- and L-resonances, the radius where 

/ e2E v y/2G F Y e X r S 



R (17) 



defines approximatively the /ir-resonance [25]. We also define here 

VL - 1 



Ar^ = tan 26*23 



/XT 



oc tan 2023?-^ (1* 



4 A first related study can been found in 
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as the resonance width, i.e the position around r Mr where the /x — r resonance effects 
take place. 

In [H], it has been shown that in the single angle approximation the neutrino- 
neutrino interaction term can be written as 

i? 4 1 

/ i ( r ) = / i 7T 2 _ fi2/r2 ( 19 ) 

where 

fi = V2G F (F*-F p R x ) (20) 

The bipolar transition starts when the synchronization ends (r = r syn ), and finish itself 
at rbip defined in [13] by the relation: 

Kru P ) * (21) 

On the left panel of Fig. ([[]), for E v = 15 MeV, the /x — r resonance occurs before the 
bipolar transition (r MT + < r bip ). Consequently, the matter eigenstates v™ and v™ 
start initially as a mix of and v T . As can be seen in Fig. ([2]), in this case the collective 
effects (represented by an arrow with a dashed line) will lead to v e m v™ where v™ is, 
at this moment, a mix of and v T . We call this type of behaviour for the parameters 
(X r ,e) type A, it is equivalent to what happens in region III in [15] . In this region, 
as in the Standard Model case, the exact place where the L-resonance occurs does not 
really modify the value of < P{y e —> u e ) > cx it- 
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Figure 1: Oscillation probabilities for a v e initially emitted as a function of the distance 
from the neutrinosphere. We show the average as dashed lines and the envelopes of 
the fast-oscillating curves as solid lines. The left, middle and right figures correspond 
respectively to type A, type B and type C behaviour with respective neutrino energy 
of 15, 30 and 35 MeV. The value of the radiative correction term is e r = 10~ 3 . 

On the right panel of Fig.flTJ, the /x — r resonance occurs right after the bipolar 
transition (r Mr — > r bip ). As can be seen in Fig. ([2]), in this case the collective effects 
(represented by an arrow with a full line) will imply that v™ w v T after the transition. 
One would have a perfect bipolar transition from v™ to z/j 71 , if the /x — r resonance 
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were well separated from the collective effects. The type C behaviour is equivalent to 
what happens for parameters (X r ,e) in region IVa in jl5]. When the fi — r resonance 
happens, it induces an exchange between the P(u e —> v^) and P(v e —> v T ) probabilities, 
which corresponds to an adiabatic evolution of z/J™ and v™, respectively (~ v e and 
~ Up before the resonance). For this reason, at the L-resonance, P(^ e —> z/ e ) strongly 
increases while P(z/ e — > u^) and P(z/ e — > u T ) decrease. 

On the middle panel of Fig.([T|), the \i — r resonance occurs during the bipolar 
transition. We call this behaviour type B and it is a transition behaviour between A 
and C. Concerning the position of the bipolar transition w.r.t. the \i — r resonance, 
we have the following inequalities: r Mr — < ru P < r^ T + ^f-. Depending on the 
exact position where the \l — t resonance happens, the bipolar transition will exchange 
v™ ~ v e mostly either with z/] 71 « v T if ru v is closer to r Mr — or with 1/3 1 fa if 
rbip is closer to r MT + ^y 1 . Such "interference" between the bipolar transition and the 
\i — r resonance is represented on Fig. ([2]) by a wiggled arrow going from v™ to a mix 
of z/j™ and v™. 

Note that the type of behaviour that neutrinos undergo depends on the parame- 
ter couple value (A r ,e) but also on the value of the initial neutrino luminosity that 
determines the strength of the collective effects. Consequently, for different values of 
parameters X r ,s, etc... the energies and/or the place at whom those behaviours occur 
will be modified. 




Figure 2: Level crossing scheme of neutrino conversion for the inverted hierarchy in a 
medium with a large V^ T with ^-mixing in the first octant. 



4 Effects of the value of e 

In this section we study the influence of the value of e, all other parameters being fixed 
to their respective values used in the previous section. In a previous paper [33] we 
showed that, according to certain value of SUSY parameters, the radiative correction 
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terms could yield values of e typically contained between ~ 5.4 x 10 -5 i.e the SM value 
and ~ 2 x 10~ 2 . Obviously the parameters space has not been fully scanned and larger 
value of e might possibly be found. Modifying this value will of course shift the position 
(see Eq JT7|) where the \i — r resonance occurs. Consequently, the energy values for a 
v e to have a behaviour type A, B or C will be different. 

Because of the particular form of the v e flux in the inverted hierarchy due to the 
neutrino-neutrino interaction, i.e the spectral split phenomenon [T71 [18], the different 
value of e may influence the flux in very different ways. In the previous section, we 
have identified three different types of behaviour for a v e going through the supernova. 
Since the \i — r resonance obviously depends on the value of the density and on the 
value of e, we can define a general quantity N^ T equal to X r e = . In this section, 
we separate the energy spectra of v e (received on Earth) in three typical regions of N^ T 
associated to typical behaviours. Here, e is supposed to be positive and the luminosity 

is L v ,r- 

4.1 Case I 

We can define this case for z/ e s that undergo a type C behaviour only for an energy 
above E u ~ 50 MeV. This is equivalent to N^ T < 1.5 x 10 6 g.cm -3 . For our density, 
the upper bound corresponds to e ~ 5 x 10 -4 . We show in Fig.(j3]) three examples for 
e = 5.4 x 10" 5 , 10~ 4 and 5 x 10~ 4 . 




E [ MeV] distance [ in units of 100km ] 



Figure 3: Left figure: v e flux on Earth as a function of the energy for three different 
values of e. Right figure: Oscillation probabilities as a function of the distance from 
the neutrinosphere with e = 5.4 x 10~ 5 and E u = 5.7 MeV. 

In the SM case, for the displayed energies, no transition can be seen for the flux in 
Figj3] 

For e = 10 -4 , the transition cannot be seen since it is above E v ~ 50 MeV. Indeed, 
at this energy, though < P{y e — > v e ) > cx it is modified by the fi — r resonance, the 
difference between the initial v e and v x fluxes is too small to produce sizeable effects 
on the final v e flux on Earth (See Eq. f[T3"|) ) . 
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For e = 5 x 1CT 4 , a type B behaviour appears at E v > 36MeV whereas the type C 
behaviour is reached for E v ~ 50 MeV. 

In the right part of Fig. ([3]), we show the oscillating probabilities for an energy below E c 
in the SM case. E c is known as the spectral split energy and here E c ~ 8 MeV. Looking 
at the electron neutrino survival probability, we see clearly the salient features of the 
v — v interaction: the synchronization regime (for r < 60km), a bipolar transition until 
the spectral split climb-up occurs (for r ~ 100km), which makes the probability goes 
back up to ~ 1 (for r ~ 200km). For this value of N^ T and the energy considered, the 
\xt resonance position is well below the synchronization regime. Therefore, only type 
A behaviour occurs. The presence of V^ T can only be seen through the fact that the 
bipolar transition occurs for both v^, and v T neutrinos as consistent with Fig.(jH) and 
the dashed arrow A. In conclusion, in this region of A^ T , the \i — t resonance happens 
before the bipolar region, hence the u e s undergo a type A behaviour as defined above 
for E v < 36MeV. 




Figure 4: Level crossing scheme of neutrino conversion for the inverted hierarchy in a 
medium with a large V^ T with ^-mixing in the first octant and an neutrino energy 
E v < E c . 

4.2 Case II 

We can define this case for u e s that undergo a type B transition for energies between 
E v ~ E c and E v ~ 50 MeV. Such behaviour corresponds to the interval 1.5 x 10 6 
g.cm -3 < iV MT < 2 x 10 7 g.cm~ 3 . For our density, the upper bound corresponds to 
e ~ 7 x 10 -3 while the lower bound is associated to e ~ 5 x 10 -4 . We show in 
Fig.© two examples of such large e value: e = 10 -3 and 2 x 10~ 3 . We also show, for 
comparison purpose, the v e flux for the SM value of e. For e = 10~ 3 , the transition 
zone appears for 28 MeV <E V <33 MeV. 

For the case e = 2 x 10~ 3 , P{y e ~^ v e) leaves a type A behaviour for E v > 10 MeV 
and the type C behaviour is reached for E v > 30 MeV. In this case, we therefore 
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have a larger transition zone, which is due to the fact that the resonance width is 
bigger. Indeed, this resonance width being proportional to r MT , it varies like (eE\ r )s, 
and therefore gets larger when e goes from 10~ 3 to 2 xl0~ 3 . This is why the flux 
moves away from the v x initial flux[j for the E v ~ 15 MeV, i.e < P{y e —> v e ) > cx it is 
close but not equal to zero. For instance, at E u = 15 MeV, < P(v e —> v e ) > cx i t =0.15 
which implies, recalling Eq.(H3), $„ e = 0.15 x (F* - F*) + F*. For E v ~ 20 MeV, 
it seems, in a misleading way that the < V[y e — > u e ) > cx it has gone down to a zero 
value again. Actually, this is due to the fact that, for this energy, the initial fluxes of 
v e and of v x cross. Consequently, the v e flux on Earth is equal to the initial flux of v x 
for E v ~ 20 MeV. Indeed, we can verify that the type B behaviour starts anew for 20 
MeV < E v < 30 MeV. 

In the right part of Fig.(jSJ), we show the oscillating probabilities for an energy 
below E c . For this value of N^, the fir resonance position is mainly below the bipolar 
transition. Nevertheless, it leads to a small disruption of the spectral split climb-up 
since the electron neutrino survival probability does not goes up back to a value of ~ 1 
but goes to ~ 0.93 instead. < P{y e — > u e ) > exit is slightly above the one yielded in 
case I. This explains the non perfect superimposition of the flux for energy below E c . 
Note that the precise interplay between the \xr resonance and the neutrino-neutrino 
interaction in this case should be investigated more deeply. 




Figure 5: Left figure: v e flux on Earth as a function of the energy for three different 
values of e. Right figure: Oscillation probabilities as a function of the distance from 
the neutrinosphere with e = 2 x 10 -3 and E u = 5.7 MeV. We show the average as 
dashed lines and the envelopes of the fast-oscillating curves as solid lines. 



4.3 Case III 

We can define this case by the interval 2 x 10 7 g.cm~ 3 < < 6 x 10 7 g.cm -3 . For 
our density, the upper bound corresponds to e ~ 2 x 10~ 2 while the lower bound is 
e ~ 7 x 10~ 3 . We show in Fig.fjSJ) two examples for e = 7.5 x 10" 3 and 2 x 10~ 2 

5 The initial v e and v x flux can be seen in Fig.©. 
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exhibiting two different cases. Note that in the case III, Ar Mr is significantly greater 
than in the previous cases I and II because it varies as oc r MT . 

In the case III, the fi — r resonance goes out of the synchronized region for E < E c 
where E c is the energy where the split occurs. Thus it will create a disruption of the 
spectral split phenomenon. 

In the case e = 7.5 x 1CT 3 , we surprisingly see that a type A behaviour occurs for E c < 
E v < 40 MeV. It turns out that the bipolar transition for theses energies is extremely 
disrupted. This may be due to a non-linear interference between the bipolar transition 
and the width of the fi — r resonance but this should deserve further investigation. For 
E > 40 MeV, we begin a transition phase. 

In the case e = 2 x 10~ 2 , for E u > E c , r Mr + > ru P , this corresponds to the 
beginning of the transition phase which lasts till r^ T — > ru v - It is important to 
note that the more the energy increases the larger the resonance width is, so the phase 
transition begins for E v < E c . Note that the spectral split signature is vanishing for 
this value of e. 




E [ MeV] distance [ in units of 100km ] 



Figure 6: Left figure: v e flux on Earth as a function of the energy for three different 
values of e. We also show here the initial v e and v x flux emitted at the neutrinosphere. 
Right figure: Oscillation probabilities as a function of the distance from the neutri- 
nosphere with e = 7.5 x 10~ 3 and E v = 5.7 MeV. We show the average as dashed lines 
and the envelopes of the fast-oscillating curves as solid lines. 

The disruption of the spectral split phenomenon is explicitly shown in the right part 
of Fig. ([6]). The bipolar transition shows huge oscillations due to a powerful interplay 
with the [it resonance, while the spectral split climb-up makes the electron survival 
probability go up only to a value of ~ 0.8 (see right part of Fig(jH])) contrary to the 
SM case where it goes back up to ~ 1. Eventually, this leads to a higher value of 
< P[y e — > u e ) > cx it which was seen on the left part of Fig. ([6]). This phenomenon is 
represented in Fig.(j4]) by the wiggled arrow B. 

Finally, note that in this section, since we studied the influence of an increasing value 
of e, we have displayed in the three cases the oscillating probabilities with E v < E c . 
Indeed, as e increases, the resonance occurs for lower and lower energies since the 
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density remains in the three cases constant. This can be seen in Eg. (1171) . It is therefore 
interesting to look at the evolution of the probabilities at the energy around which the 
fj, — t resonance occurs. The evolution probabilities for E v > E c can be seen in Figs. 



5 The degeneracy of e with respect to other param- 
eters 

As seen in Sec. [3j the typical behaviour A, B or C are the consequence of a certain 
position of the \x — r resonance in comparison with the bipolar transition. In the 
following subsections, we study in particular one parameter that influence either the 
resonance or the v — v interaction, the other being fixed to their reference value as 
defined in Sec. [2j 



5.1 Effects of #23 and of the sign of e 

Changing the value of #23 in such way that it changes of octant, will automatically 
prevent a \i— r resonance as discussed in J5U1I1S]. Using the two flavour MSW condition, 

A 2 

cos 29 23 = V = eV2G F Y ePB , (22) 

one can see that changing the sign of e could be immediately interpreted as perfectly 
equivalent to changing of octant. However, in the case of a negative e = —\e\, looking 
at the matter interaction Hamiltonian, we have: 

H m = diag(V e ,0,-\V^ T \) = diag{V e + \V^ T \, \V^ T \,0). (23) 

In this case, the term |V^ T | becomes a matter potential for which can be seen like 
an effective presence of muon in the supernova whereas the electron matter potential 
appears slightly modified of the order of e, i.e 2% at most. 

This explains why, looking at the left figure panel in Fig. (JTj) , the top and lower part 
show opposite behaviours for and v T in the first 100km; The same thing happens for 
the right figure panel. Even when the MSW condition is not fulfilled, a matter effect 
is still present as can be seen in [5U]. This is called the matter suppression effect and 
a transition, though partial is happening. 
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Figure 7: Oscillation probabilities as a function of the distance from the neutrinosphere. 
We show the average as dashed lines and the envelopes of the fast-oscillating curves as 
solid lines. The figures correspond to the value of the couple of parameters (e = ±e 2 , 
6*23 = 40° or 50°). Si is still equal to 2 x 1CT 3 . In these figures, the neutrino energy is 
30 MeV. 

In Fig.flEl), we observe that the energy spectra of the electron neutrino flux on 
Earth, is the same for the parameter couples (e > 0, #23 = 40°) and [e < 0, #23 = 50°) 
whereas a identical shift for the transition around 30 MeV is observed in the cases 
(e < 0, #23 = 40°) and (e > 0, #23 = 50°). The fact that a similar behaviour between 
the cases (e > 0, #23 = 40°) and (e < 0, #23 = 50°) is observed is due to the fulfillment 
of the MSW condition (Eq. (l22l) ) in the two cases. The only differences can be seen on 
the evolution with the distance of the and v T oscillation probabilities, as said above 
because in the first/second case, the neutrinos see respectively an effective presence 
of r/fj, leptons. In the opposite cases, no MSW resonance can happen but a matter 
suppression effect occurs. To explain this shift in energy, one can think the matter 
suppression as a transition between and v T neutrinos but less effective as the MSW 
resonance. Since in all cases, all the parameters except the sign of e and the octant of 
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6*23, are identical, to maximize the effect of this matter suppression, one needs to move 
it away from the collective effects region, this can be done by increasing the neutrino 
energy. Note that the matter suppression can interfere with the collective effects as 
represented as a wiggled arrow in Fig.(JH]). 




„K 1 1 1 L_ 

10 20 30 40 

E [ MeV] 



Figure 8: v e flux on Earth as a function of the energy for different values of the 
parameter couple ( e, 623). e 2 is still equal to 2 x 10 -3 . 




Figure 9: Level crossing scheme of neutrino conversion for the inverted hierarchy in a 
medium with a large V^ T with ^23-niixing in the second octant. 

5.2 Effects of the luminosity 

Effects of luminosity should not be neglected. As the cooling phase passes, the neutrino 
fluxes emitted at their respective neutrinosphere, decrease in intensity. Such diminution 
can be approximated by a decreasing exponential as seen in Eq. ffT2"j) . Consequently, the 
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strength of the v — v interaction shrinks, implying a slow diminution of the collective 
effect predominating zone: 

r bip ~ $ ~ exp f-^-J (24) 

As the bipolar radius decreases, r MT being fixed, the electron neutrino probability be- 
haviour tends to be like the type C behaviour. This can be seen in Fig.f flOl) where the 
fluxes have been normalized to the reference value in order to appreciate only the \i — r 
j v — v interplay effects. Knowing that r^ T ~ (A r e)3, multiplying the luminosity by 

3 

a factor a would in principle be equivalent to multiply either A r or e by a factor a~ 
since 

^ ~ (X r e)hf (25) 

Therefore, increasing the luminosity or shifting inwards the position of the \i — t reso- 
nance is rather equivalent until the neutrino density is too weak to maintain collective 
effects. 

Looking at the left part of Fig. ffTO]) , one sees the energy spectra of v e varying as 
predicted above. When L v = 2 x 10 51 erg.s _1 the v e present a similar shape with the 
reference case, except that the transition zone starts for E v ~ 35 MeV. Indeed, a higher 
luminosity is equivalent to a \x — r resonance position deeper inside the supernova. 
Therefore, neutrinos need a higher energy value to begin to see the type B behaviour. 
Symmetrically, a lower value of the initial neutrino flux will make the type B behaviour 
appears for lower energies, i.e E v ~ 23 MeV in this case. Note that because of a low 
value of the neutrino flux, the collective effects start to weaken which leads to some 
oscillations in the flux for energies E v < E c . Equivalently, one can look at the right 
part of Fig.f llOp where we show the electron neutrino survival probability for the three 
different values of the luminosity: L u ^ r , L Vj i and L Ut2 - The value of < P(u e —> u e ) > cx it 
is anew highly dependent of the relative position of r Atr w.r.t rup- This figure shows 
explicitly the increase of r^p with the luminosity. 
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Figure 10: Left figure: luminosity-normalized v e flux on Earth as a function of the 
energy for three different luminosities. Right figure: Electron neutrino survival proba- 
bilities as a function of the distance from the neutrinosphere for three different lumi- 
nosities with E v = 30 MeV. We show the average as dashed lines and the envelopes of 
the fast-oscillating curves as solid lines. 



5.3 Effects of the density 

Depending on the mass of the progenitor star, the explosion releases a variable energy 
which, under the form of a shock-wave, can lead to important variation on the matter 
density, whose precise profile remains partially unknown. In our analytical profile, it 
means that the value of A r will decrease with time. The effect of the density variation 

through the parameter A r is crystal clear. Using Eq. flTT]) . one can see that the position 

i 

of the n — t resonance varies like Ar . Therefore, changing the density or the value 
of the radiative correction term e is almost perfectly equivalent as long as the density 
does not reach very high value such as in the accretion phase. This can be seen in 
the left part of Fig. ffTT]) where the density is Ai = 2A r and A2 = A r /2 are respectively 
equivalent to the reference case but with E\ = 2 x 10 -3 and 62 = 5 x 10~ 4 . The 
incertitude concerning the density profile plays an important role as gives rise to 
three well distinct behaviours corresponding to the cases I, II and III. Equivalently, 
one can look at the right part of Fig. ffTT]) . where we show the electron neutrino survival 
probability for the three different values of the density: A r , Ai and A2; We can notice 
that rbip is the same for the three curves, which is consistent with the fact that the 
luminosity has not been modified. The bipolar transition is a bit disrupted in the case 
where the matter is more dense. This is due to the fact that the V^ r term is larger, 
which implies a more powerful interplay with the fi — r resonance. Moreover, we see 
a light shift of the L-resonance position towards the outside of the supernova with an 
increasing value of the density parameter A. 
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Figure 11: Left figure: v e flux on Earth as a function of the energy for three different 
densities. Right figure: Electron neutrino survival probabilities as a function of the 
distance from the neutrinosphere for three different densities with E v = 30 MeV. We 
show the average as dashed lines and the envelopes of the fast-oscillating curves as 
solid lines. 



6 Detection of the signal on Earth 

Since the explosion of SN1987A in the Large Magelanic Cloud, it has been proved 
that neutrinos from supernova can be detected on Earth. Such astrophysical event 
can yield a tremendous amount of information through the neutrino fluxes received on 
Earth. With the currently running neutrino observatories, a high-statistics signal for a 
galactic SN explosion is expected in the electron anti-neutrino channel. For instance, 
Super-Kamiokande could detect through the reaction v e + p — > n + e + about ~ 10 4 
events. Nevertheless, larger detectors could yield an even better statistics and therefore 
could be sensitive to more discreet features of the explosion. Also, they logically could 
detect a more distant explosion with the same current number of event statistics. 

Three main types of multi-purpose detectors are proposed to enhance the quality 
of information about the neutrino properties and the explosion mechanism that we will 
receive from the supernova neutrino fluxes. The proposed future neutrino observatories 
offered complementary detection techniques: a megaton water Cherenkov detector such 
as MEMPHYS or Hyper-Kamiokande, a liquid scintillator like LENA or a liquid argon 
Time Projection Chamber (TPC) like GLACIER. See e.g. jl6] for a detailed review 
on those large-scale future detectors. The interest of such detectors, for a supernova 
neutrino signal, has been showed e.g. in jlQ] for the water-Cherenkov type or in [37] 
for a liquid Argon type. We focus on a detector of the latter kind since it has the 
possibility to detect electron neutrino by charged- current, through the reaction 

v e + 40 Ar ^ 40 K* + e~ (26) 

therefore yielding a very clear signal complementary of the inverse /3 of the water 
Cherenkov detectors. 
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To calculate the number of events received in such liquid Argon detector (Icarus- 
like), we use the cross section given in [37] for the charged- current reaction (126]) . For 
the detection of the \i — r resonance effect in a liquid Argon type experiment, we 
calculate the usual differential number of neutrinos detected at a distance D from the 
supernova: 



where N ta rgets is the number of targets in the detector, D is the distance from the 
supernova, E v is the neutrino energy, t is the detector time, d 2 N u / ' dE v dt the differential 
flux of neutrinos arriving on Earth, for a given detector time and a given neutrino 
energy, and <j(E u ) is the cross section of the reaction process considered. For a galactic 
explosion, we suppose that the supernova explosion will happen at the distance D = 
lOkpc. Our calculations do not include Earth matter effects and the Poisson error 
due to the finite number of detected events. 

For the signal really obtained in the detector, note that the energy resolution 
reached in the detector is good enough to appreciate precisely the fi — r resonance 
effect. Indeed, the energy of the neutrino will not be directly obtained. It will be 
accessible via the electron energy of the reaction of Eq. (I26j) . 

In a realistic supernova environment, after the post-bounce, both the density profile 
and the luminosity are decreasing. As seen in the previous section, this affects the 
interplay between the \x— r resonance and the v—v interaction. If the density decreases, 



the fj, — t resonance tends to get closer from the nascent iron core with r^ T ~ . 



When the luminosity decreases, the bipolar region shrinks as r bip ~ L$. Those two 
phenomena are therefore opposite and tend to provoke a status quo concerning the 
type of behaviour of the electron neutrino flux. The way the v e flux energy distribution 
is modified with time depends on the relative variation speed of the parameters A r and 

Of course, the absolute value of the flux will undergo a global exponential decrease 
with time. In the literature tackling the question of dynamical density profile induced 
by shock-waves, the difference of density during a time variation of about 5s, tends to 
be of a factor ~ 2, for the part of the density which can be well approximated by a 
power-law. See e.g. [H] for the density profile as a function of time. In Figs. (|T2|) and 
(fl3|) . we show the differential number of events received in a 70kton liquid argon type 
detector and the associated flux for a post-bounce time of 3 and 8 s. We study three 
distinct values for the radiative correction term e: the SM value Ssm = 5.4 x 1CT 5 whose 
N^r will stay in region I for these two times, an intermediate value e r = 1CT 3 which 
corresponds for the v e flux to a case II behaviour until 8s, and the typical maximum 
value obtained in the SUSY framework, i.e £4 = 2 x 10~ 2 which corresponds for the v e 
flux to case III behaviour during most of the cooling phase. 

6 Because they depend upon the the position of the supernova with respect to the detector when 
the event occurs, we do not consider them in the neutrino signal in the detector. A study of their 
effects is given for instance in [T7] . 




(27) 
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Figure 12: In this case, the time after post-bounce is supposed to be t=3s, the initial 
luminosity taken is L u (t = 3s) = 4.24x 10 51 erg.s -1 , the density parameter is A r = 1.14 
x 10 6 and 6^3 = 40°. Left figure: differential number of events on Earth in a Liquid 
Argon detector of 70 kton size as a function of the energy for three different e. Right 
figure: corresponding v e flux on Earth as a function of the energy for three different e. 

For t = 3s, the neutrino fluxes with e r or e 4 entering the detector show important 
differences with the SM case. While for e r the difference with esm appears for an 
energy of ~ 32MeV as a sizeable decrease of flux, the difference for £4 appears at a 
lower energy ~ 22MeV. Looking at the differential number of events, we see that the 
usual spectral split feature becomes almost invisible except in the case of £4 where a 
very small bump remains. On the contrary the effects of the \i — r resonance at higher 
energies are much more visible. In the case of e r , a new spectral split feature appears. 
A type B behaviour starts at ~ 32 MeV and leads to a type C behaviour from ~ 38 
MeV. For £4, the flux undergoes a diminution from ~ 22MeV but this split feature 
is not as salient as the previous one. Nevertheless, we see that in the £4 case, the 
number of events is globally smaller than the e r case. With respect to the SM case, 
the difference in the number of events corresponds in average to a factor ~2 from ~ 22 
MeV. 
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Figure 13: In this case, the time after post-bounce is supposed to be t=8s, initial 
luminosity taken is L u (t=8s) = 10 51 erg.s -1 , the density parameter is A r = 5.7 x 10 5 
and 823 = 40°. Left figure: differential number of events on Earth in a Liquid Argon 
detector of 70 kton size as a function of the energy for three different e. Right figure: 
corresponding v e flux on Earth as a function of the energy for three different e. 

For t = 8s, the neutrino fluxes with e r or e 4 display a modification w.r.t the SM case 
similarly to time t = 3s. Concerning the flux with e r , it leaves the type A behaviour for 
a lower energy equals to ~ 28 MeV. this can be understood recalling Eq. (T25|) . Indeed, 
since A r has been multiplied by a factor 0.5 and the luminosity has decreased by a 
factor ~ 4, the relative position of r Atr over r^p has increased by a factor ~ 43 2~3 > 1. 
Consequently, a lower energy is required for the type B behaviour to start. Concerning 
the flux with £4, the relative difference w.r.t. the flux in the SM case is more important 
at low energies. Nevertheless, the luminosity being exponentially decreasing with time 
and the low value of the cross-section at low energies will prevent to see the effects of 
£4 on the number of events for those energies. Finally, we can see on the number of 
events figure, that the e r curve presents globally a more pronounced difference w.r.t to 
the SM case than the £ 4 curve. 
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Figure 14: In this case, the time after post-bounce is supposed to be t=8s, initial 
luminosity taken is L u (t=8s) = 10 51 erg.s -1 , the density parameter is A r = 5.7 x 10 5 
and # 2 3 = 40°. In these figures the hierarchy in temperature is (E Ve ) = 12 MeV, 
(E Pe ) = 15 MeV, (E Vx ) = 18 MeV. Left figure: differential number of events on Earth 
in a Liquid Argon detector of 70 kton size as a function of the energy for three different 
e. Right figure: corresponding v e flux on Earth as a function of the energy for three 
different e. 

In Fig.(ll4p. we take other values for the energy hierarchy in agreement with the 
numerical simulations performed in (45]: i.e. (E Ue ) = 12 MeV, (E Pe ) = 15 MeV, 
(E Ux ) = 18 MeV respectively. Concerning the flux with e r , the difference w.r.t the SM 
case appears from ~ 23 MeV. For £4, the spectral split feature has vanished and the 
flux presents a much higher value w.r.t the SM case around 10 MeV. Concerning the 
number of events, both e r and £4 curves move away from the SM one from around 20 
MeV. We therefore conclude that even with this hierarchy of temperatures, sizeable 
effects are present on the differential number of events received in a Liquid Argon 
detector. 

Though we did not look at the time evolution v e flux for given energies during the 
whole cooling phase, it seems that the global pattern of the fluxes will not be much 
modified with time, contrary to what may happen for the electron anti-neutrino flux 
|34j . In this case, to distinguish w.r.t. the SM case would require the identification of 
a particular energy distribution of the flux. 

Considering the v e flux, for the e r case, the spectral split feature seems peculiar 
enough to be recognized: at 40 MeV, the ratio between the SM flux and the e r flux is 
more or less constant. One could also look at the time evolution of the flux for an energy 
of 30 MeV and 40 MeV. In the SM case, the flux at 40 MeV will be systematically 
higher than the one at 30 MeV. In the case of a large e, the flux at 40 MeV will be 
systematically lower. For the £4 case, the low energy bump may be visible for early 
time because of an important global luminosity. Therefore, the presence of a spectral 
split at high energy in the differential number of events may be a signature of a large 
radiative correction term and therefore a beyond standard physics. 
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Using the signal in the electron anti-neutrino channel with a Cherenkov detector, 
one would be able to obtain precise information on the density profile as a function of 
time. In addition, with enough precision on the neutrino parameters, like the octant 
of #23; as well as Supernova parameters partly given by the other channels (electron 
anti-neutrino via inverse /3 reaction, all flavour flux via neutral current), this could 
ideally lead to new constraints for the SUSY parameters via the analytical calculation 
of the radiative correction as in [331 [23]. Equivalently, the absence of such characteristic 
imprints would also lead to constraints if no other effects like turbulence modify the v e 
flux existing inside the supernova as will be discussed in the next section. 

Finally, note that the impact of the interplay between V^ T and the collective effects 
has been observed on the Diffuse Supernova Neutrino Background (DSNB) [53]. Such 
interplay tends to provoke important effects even in the Standard Model thanks to the 
redshift. 

7 Discussions 

In this section we discuss on how much these results can be seen as realistic. 

We have performed a calculation of the neutrino propagation in a core-collapse su- 
pernova including the neutrino-neutrino interaction and a large mu-tau index of refrac- 
tion coming from SUSY radiative corrections. Such combination in certain conditions 
of matter density, energy, luminosity and oscillation parameters implies a resonance, 
associated to the loop correction potential. This occurs in the same region as the bipo- 
lar resonance. Consequently, the and v T oscillation probabilities are modified and 
the low-resonance becomes much more adiabatic leading to an important inversion of 
fluxes, visible on the electron neutrino flux on Earth as a function of energy. Such 
characteristic imprints of a large V^ T , i.e of a BSM physics, would be also sizeable on 
the number of events in a Liquid Argon detector type. 

In a supernova environment, such interplay has been first studied for a given energy 
in [30] . In order to obtain a fi — r resonance with a radiative correction term in the SM 
framework, the authors considered a very high density profile. This way, such resonance 
happens well outside the bipolar transition zone and a clear, easily understandable, 
behaviour is exhibited. Unfortunately, it has been shown in [21], for the multi-angle 
and, therefore, general case, that such high density prevents collective effect from 
taking place. Consequently, we have to check if the density we use is low enough for 
the multiangle matter effect to be neglected. In [15], they found that the condition: 



must be fulfilled. Using Eq.(ll) of [H], we find in our case, fiQ ~ 2.6 x 10 5 km" 1 . 
Therefore, having Y e = 0.5 and assuming r syn ~ 100 km the condition becomes : 
A r < 5.2 x 10 6 km" 1 . Since A r = 5.7 x 10 5 km -1 , we can indeed safely neglect the 
matter effects. 




syn 
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The advantage of having a large V^ T potential due to a large SUSY radiative correc- 
tion is significant for multiple reason. In this case, a very high density is not required 
for the \i — t resonance to be outside the zone where collective effects predominate. 
Therefore, we naturally avoid the problem of the high density effect cited above. 

Concerning our density profile, we use an analytical power law which does not take 
into account shock-waves effects. Nevertheless we believe such dynamical effects are 
not determinant in our case. Indeed, we checked [50] that the shock-wave effects such as 
depicted in [33] do not affect the v e In inverted hierarchy, z/ e s do not encounter the 

H-resonance, and therefore do not undergo the effect of adiabaticity breaking and/or 
multiple resonance which tends to smear out the neutrino flux distribution. Moreover, 
looking at the v e flux during the cooling phase not too early, allows the shock wave 
to affect the density profile well after 1000 km, i.e well after that the \i — t resonance 
and the v — v interaction occur. Similarly, the L resonance is not supposed to be 
significantly affected by the shock- wave effects [4"T] . 

Concerning the collective effects, even if the matter density is not so high, multi- 
angle should be prevented of triggering decoherence because of the asymmetry between 
the v e and v e flux cause by deleptonization [33]- Consequently, single-angle behaviour 
may well be typical for realistic SN conditions as long as this flux asymmetry remains. 

Note that we have checked that the influence of S is negligible in this particular case 
where we focus on the electron neutrino signal in the inverted hierarchy. Moreover, we 
took equal luminosities for and v T flux [27]. As long as it is non zero, the influence of 
#13 on the v e flux in inverted hierarchy is very small, possibly below the experimental 
uncertainties. 

In our paper, we also made the assumption of a constant electron fraction. This 
approximation is justified when we restrict ourselves to a maximum value of ~ 2 x 10~ 2 
for radiative term e. In this case, no non-standard interaction (NSI) resonance will 
occur as shown in [35]. |^| 

Similar interplay should be in principle present in the anti-neutrino channel but as 
pointed out in [33], the presence of shock- waves effects tend to render the H-resonance 
non adiabatic and eventually create a multiple resonance effect which can smear out 
the electron anti-neutrino signal. The addition of such important radiative correc- 
tion matter potential should disturb the signal but full numerical calculation as first 
performed in should be done to witness the impact on the electron anti-neutrino flux. 

Concerning the normal hierarchy, the effects of the \i — r resonance tend to be 
much smaller as no bipolar transition occurs. Nevertheless, the case where spectral 
splits take place in this hierarchy [511 [52] should be investigated in the presence of a 
large s. 

Finally, we discuss the most important effect that could prevent the salient modi- 

7 The results found in [35] show explicitly that the presence of the Shockwave does not prevent the 
consequences on the electron neutrino flux of the interplay between V^ T and the collective effects. 

8 Note that our value of e doesn't match exactly to the value of the NSI diagonal term e TT given 
in gSJ. 
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fication of the v e flux induced by the \x — r resonance and v — v interplay. A realistic 
density profile should present possible stochastic matter fluctuations due to hydrody- 
namical instabilities created in the wake of a shock front. This has been studied in 
several works [521 EI] but the precise size of such fluctuations remains an important 
open issue. In principle, the \i — t resonance should not be affected. The question 
whether the L-resonance will be affected or not, depends on the fluctuation scale i.e if 
the damping effects are significant on the flavour transition pattern. Associated with 
this issue, taking into account a non-spherical symmetry for the neutrino emission [55] 
as well with matter density profile [JT], may be relevant, for example, in the presence 
of hydrodynamic turbulence in the background of ordinary matter. The framework 
considered in this paper, defined in Sec. [21 allows us to think that our results should 
contain the main features of a realistic signal. 



8 Conclusion 



In this paper, we worked in the SUSY framework. If new phenomena compatible 
with this framework are discovered on the supernova neutrino fluxes, one of the major 
challenges will be to find out the underlying model and to measure its parameters as 
precisely as possible. This can be done for example by using codes in which we consider 
the phenomenological constraints [56] . It will be necessary to investigate the precisiorj^] 
with which SUSY model parameters [33] can be derived from these measurements. 
Finally, if Supersymmetry is realised in Nature, it should be possible to constrain the 
SUSY parameters thanks to the discovery of sparticles at TeV colliders. For a galactic 
explosion, using a fully detailed numerical simulation, supernova neutrino fluxes could 
ideally lead to constraints on the SUSY parameters in the inverted hierarchy case. In 
the future, the results presented in this paper should be applied to a larger variety of 
models and confronted to some other data from the Large Hadron Collider (LHC) and 
the International Linear Collider (ILC). Finally, let us say that the /i — r resonance 
effect on the energy distribution of the electron neutrino flux could impact on the 
number of events for relic supernova neutrinos as in [35]. In this case, it could give an 
even lower total red-shift-integrated number of events than in the SM case which can be 
consistent with observations. Moreover, a large e could also impact the heavy elements 
nucleosynthesis in supernova, and more precisely on the electron fraction through the 
r-process. Supernovae can definitively be seen more than ever, as an open window 
towards beyond Standard Model physics. 
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